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Abstract 

We propose an entropic geometrical model of psycho-physical crowd dynamics (with dis- 
sipative crowd kinematics), using Feynman action-amplitude formalism that operates on 
three synergetic levels: macro, meso and micro. The intent is to explain the dynamics 
of crowds simultaneously and consistently across these three levels, in order to character- 
ize their geometrical properties particularly with respect to behavior regimes and the state 
changes between them. Its most natural statistical descriptor is crowd entropy S that sat- 
isfies the Prigogine's extended second law of thermodynamics, dtS > (for any nonisolated 
multi-component system). Qualitative similarities and superpositions between individual 
and crowd configuration manifolds motivate our claim that goal-directed crowd movement 
operates under entropy conservation, dtS — 0, while natural crowd dynamics operates under 
(monotonically) increasing entropy function, dtS > 0. Between these two distinct topological 
phases lies a phase transition with a chaotic inter-phase. Both inertial crowd dynamics and 
its dissipative kinematics represent diffusion processes on the crowd manifold governed by 
the Ricci how, with the associated Perelman entropy-action. 

Keywords: Crowd psycho-physical dynamics, action-amplitude formalism, crowd mani- 
fold, Ricci flow, Perelman entropy, topological phase transition. 
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1 Introduction 



Recall that the term cognitioT^ is used in several loosely related ways to refer to a faculty 
for the human-like processing of information, applying knowledge and changing preferences. 
In psychology, cognition refers to an information processing view of an individual psycho- 
logical functions (see ^651 [771 [Ml |3l [2 ) ■ More generally, cognitive processes can be natural 
and artificial, conscious and not conscious; therefore, they are analyzed from different per- 
spectives and in different contexts, e.g., anesthesia, neurology, psychology, philosophy, logic 
(both Aristotelian and mathematical), systemics, computer science, artificial intelligence 
(AI) and computational intelligence (CI). Both in psychology and in AI/CI, cognition refers 
to the mental functions, mental processes and states of intelligent entities (humans, human 
organizations, highly autonomous robots), with a particular focus toward the study of com- 
prehension, inferencing, decision-making, planning and learning (see, e.g. [lOj). The recently 
developed Scholarpedia, the free peer reviewed web encyclopedia of computational ncuro- 
science is largely based on cognitive neuroscience (see, e.g. [7S]). The concept of cognition is 
closely related to such abstract concepts as mind, reasoning, perception, intelligence, learn- 
ing, and many others that describe numerous capabilities of the human mind and expected 
properties of AI/CI (see [SS1I31] and references therein). 

Yet disembodied cognition is a myth, albeit one that has had profound influence in West- 
ern science since Rene Descartes and others gave it credence during the Scientific Revolution. 
In fact, the mind-body separation had much more to do with explanation of method than 
with explanation of the mind and cognition, yet it is with respect to the latter that its 
impact is most widely felt. We find it to be an unsustainable assumption in the realm of 
crowd behavior. Mental intention is (almost immediately) followed by a physical action, that 
is, a human or animal movement [78] . In animals, this physical action would be jumping, 
running, flying, swimming, biting or grabbing. In humans, it can be talking, walking, driv- 
ing, or shooting, etc. Mathematical description of human/animal movement in terms of the 
corresponding neuro-musculo-skeletal equations of motion, for the purpose of prediction and 
control, is formulated within the realm of biodynamics (see [4T1 [53l l42l l43l l44l l45l [46l I47j ) . 

The crowd (or, collective) psycho-physical behavior is clearly formed by some kind of su- 
perposition, contagion, emergence, or convergence from the individual agents' behavior. Le 
Bon's 1895 contagion theory, presented in "The Crowd: A Study of the Popular Mind" influ- 
enced many 20th century figures. Sigmund Freud criticized Le Bon's concept of "collective 
soul," asserting that crowds do not have a soul of their own. The main idea of Freudian crowd 
behavior theory was that people who were in a crowd acted differently towards people than 
those who were thinking individually: the minds of the group would merge together to form 
a collective way of thinking. This idea was further developed in Jungian famous "collective 
unconscious" [ST]. The term "collective behavior" [7] refers to social processes and events 
which do not reflect existing social structure (laws, conventions, and institutions), but which 
emerge in a "spontaneous" way. Collective behavior might also be defined as action which 
is neither conforming (in which actors follow prevailing norms) nor deviant (in which actors 
violate those norms) . According to the emergence theory [82 , crowds begin as collectivities 
composed of people with mixed interests and motives; especially in the case of less stable 

^Latin: "cognoscere — to know" 
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crowds (expressive, acting and protest crowds) norms may be vague and changing; people in 
crowds make their own rules as they go along. According to currently popular convergence 
theory, crowd behavior is not a product of the crowd itself, but is carried into the crowd by 
particular individuals, thus crowds amount to a convergence of like-minded individuals. 

We propose that the contagion and convergence theories may be unified by acknowledging 
that both factors may coexist, even within a single scenario: we propose to refer to this 
third approach as behavioral composition. It represents a substantial philosophical shift 
from traditional analytical approaches, which have assumed either reduction of a whole into 
parts or the emergence of the whole from the parts. In particular, both contagion and 
convergence are related to social entropy, which is the natural decay of structure (such as 
law, organization, and convention) in a social system [15j . Thus, social entropy provides an 
entry point into realizing a behavioral-compositional theory of crowd dynamics. 

Thus, while all mentioned psycho-social theories of crowd behavior are explanatory only, 
in this paper we attempt to formulate a geometrically predictive model-theory of crowd 
psycho-physical behavior. 

We propose the entropy formulation of crowd dynamics as a three-step process involving 
individual psycho-physical dynamics and collective psycho-physical dynamics. 



2 Generic three— step crowd behavioral dynamics 

In this section we propose a generic crowd psycho-physical dynamics as a three-step process 
based on a general partition function formalism. Note that the number of variables Xi in 
the standard partition function from statistical mechanics (see equation (jSSp in Appendix) 
need not be countable, in which case the set of coordinates {a;*} becomes a field cj) = <^(a;), 
so the sum is to be replaced by the Euclidean path integral (that is a Wick-rotated Feynman 
transition amplitude in imaginary time, see subsection 13. 4p . as 

Zi^) = J V[4>] exp [-77(0)] . 

More generally, in quantum field theory, instead of the field Hamiltonian H{4)) we have the 
action S{4)) of the theory. Both Euclidean path integral, 

Z{<j)) = j 'D[<j)\ exp [— S'((/')] , real path integral in imaginary time (1) 

and Lorentzian one, 

Z{(j)) = j exp [^^(i?!))] , complex path integral in real time (2) 

-r epresent quantum field theory (QFT) partition functions. We will give formal definitions 
of the above path integrals (i.e., general partition functions) in section 3. For the moment, 
we only remark that the Lorentzian path integral ^ represents a QFT generalization of the 
(nonlinear) Schrodinger equation, while the Euclidean path integral ([T]) in the (rectified) real 
time represents a statistical field theory (SFT) generalization of the Fokker-Planck equation. 
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Now, following the framework of the Extended Second Law of Thermodynamics (see 

Appendix), diS > 0, for entropy S in any complex system described by its partition function, 
we formulate a generic crowd psycho-physical dynamics, based on above partition functions, 
as the following three-step process: 

1. Individual psycho physical dynamics (I'D) is a transition process from an entropy- 
growing "loading" phase of mental preparation, to the entropy-conserving "hitting/shooting" 
phase of physical action. Formally, 1T> is given by the phase-transition map: 

"LOADING" ■.dtS>0 "HITTING" : dt S=0 



IV: MENTAL PREPARATION =^ PHYSICAL ACTION (3) 
defined by the individual (chaotic) phase-transition amplitude 



/PHYS. ACTION 



CHAOS 



dtS>0 \ f r, r*i 

MENTAL PREP.) := / ©[^j e*'^™!*! 
/id J 



where the right-hand-side is the Lorentzian path-integral (or complex path-integral in real 
time, see Appendix), with the individual psycho-physical action 

Slum = I '^^ LiY^i^dt, 



where Lid ["I^] is tlu^ psycho-physical Lagrangian, consisting of mental cognitive potential and 
physical kinetic energy. 

2. Aggregate psycho-physical dynamics {AD) represents the behavioral composition- 
transition map: 



"LOADING" : (9t S>0 "HITTING" : dt S=0 



AV : ^ MENTAL PREPARATION =^ ^ PHYSICAL ACTION^ (4) 

ieAD ieAD 

where the (weighted) aggregate sum is taken over all individual agents, assuming equipar- 
tition of the total psycho-physical energy. It is defined by the aggregate (chaotic) phase- 
transition amplitude 



/ dtS=o dtS>o \ r 

^PHYS. ACTION CHAOS MENTAL PREP.^ := J V[^] e 



-Sad[*] 



with the Euclidean path-integral in real time, that is the SFT-partition function, based on 
the aggregate psycho-physical action 

5ad[$]=/ LAD[<^]dt, with Lad[$] = ^ iro[^]- 

Jtini ieAU 

3. Crowd psycho-physical dynamics {CV) represents the cumulative transition map: 

"LOADING" ■.dtS>0 "HITTING" : dt S=0 

CV : ^ MENTAL PREPARATION ^ PHYSICAL ACTION^ (5) 

ieCD ieCD 
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where the (weighted) cumulative sum is taken over all individual agents, assuming equiparti- 
tion of the total psycho-physical energy. It is defined by the crowd (chaotic) phase-transition 
amplitude 



( 



dtS=0 dtS>0 

PHYS. ACTION CHAOS MENTAL PREP. 



) 



CD 



/ 



with the general Lorentzian path-integral, that is, the QFT-partition function), based on 
the crowd psycho-physical action 




ieCD fc= # of ADs in CD 



All three entropic phase-transition maps, TV, AV and CD, are spatio-temporal bio- 
dynamic cognition systems, evolving within their respective configuration manifolds (i.e., 
sets of their respective degrees-of- freedom with equipartition of energy) , according to bipha- 
sic action- functional formalisms with psycho-physical Lagrangian functions Lid, -^ad and 
LcD, each consisting of: 

1. Cognitive mental potential (which is a mental preparation for the physical action), and 

2. Physical kinetic energy (which describes the physical action itself). 

To develop I'D, AD and CD formalisms, we extend into a physical (or, more precisely, 
biodynamic) crowd domain a purely-mental individual Life-Space Foam (LSF) framework 
for motivational cognition [52 , based on the quantum-probability concept o 

■^The quantum probability concept is based on the following physical facts |56l [57] 

1. The time-dependent Schrodinger equation represents a complex-valued generalization of the real-valued 
Fokker-Planck equation for describing the spatio-temporal probability density function for the system 
exhibiting continuous-time Markov stochastic process. 

2. The Feynman path integral (including integration over continuous spectrum and summation over discrete 
spectrum) is a generalization of the time-dependent Schrodinger equation, including both continuous-time 
and discrete-time Markov stochastic processes. 

3. Both Schrodinger equation and path integral give 'physical description' of any system they are modelling in 
terms of its physical energy, instead of an abstract probabilistic description of the Fokker-Planck equation. 

Therefore, the Feynman path integral, as a generalization of the (nonlinear) time-dependent Schrodinger equation, 
gives a unique physical description for the general Markov stochastic process, in terms of the physically based 
generalized probability density functions, valid both for continuous-time and discrete-time Markov systems. Its 
basic consequence is this: a different way for calculating probabilities. The difference is rooted in the fact that 
sum of squares is different from tie square of sums, as is explained in the following text. Namely, in Dirac- 
Feynman quantum formalism, each possible route from the initial system state A to the final system state B is 
called a history. This history comprises any kind of a route, ranging from continuous and smooth deterministic 
(mechanical-like) paths to completely discontinues and random Markov chains (see, e.g., [22]). Each history 
(labelled by index i) is quantitatively described by a complex number. 

In this way, the overall probability of the system's transition from some initial state A to some final state B is 
given not by adding up the probabilities for each history-route, but by 'head-to-tail' adding up the sequence of 
amplitudes making-up each route first (i.e., performing the sum-over-histories) - to get the total amplitude as a 
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The psycho-physical approach to XT), AD and CD is based on entropic motor control 
[391 140] , which deals with neuro-physiological feedback information and environmental un- 
certainty. The probabilistic nature of human motor action can be characterized by entropies 
at the level of the organism, task, and environment. Systematic changes in motor adap- 
tation are characterized as task-organism and environment-organism tradeoffs in entropy. 
Such compensatory adaptations lead to a view of goal-directed motor control as the product 
of an underlying conservation of entropy across the task-organism-environment system. In 
particular, an experiment conducted in [40j examined the changes in entropy of the coordi- 
nation of isometric force output under different levels of task demands and feedback from 
the environment. The goal of the study was to examine the hypothesis that human motor 
adaptation can be characterized as a process of entropy conservation that is reflected in the 
compensation of entropy between the task, organism motor output, and environment. Infor- 
mation entropy of the coordination dynamics relative phase of the motor output was made 
conditional on the idealized situation of human movement, for which the goal was always 
achieved. Conditional entropy of the motor output decreased as the error tolerance and 
feedback frequency were decreased. Thus, as the likelihood of meeting the task demands was 
decreased increased task entropy and/or the amount of information from the environment 
is reduced increased environmental entropy, the subjects of this experiment employed fewer 
coordination patterns in the force output to achieve the goal. The conservation of entropy 
supports the view that context dependent adaptations in human goal-directed action are 
guided fundamentally by natural law and provides a novel means of examining human motor 
behavior. This is fundamentally related to the Heisenberg uncertainty principle [57j and 
further supports the argument for the primacy of a probabilistic approach toward the study 
of biodynamic cognition systems H 

'resultant vector', and then squaring the total amplitude to get the overall transition probability. 

Here we emphasize that the domain of validity of the 'quantum' is not restricted to the microscopic world [83] . 
There are macroscopic features of classically behaving systems, which cannot be explained without recourse to 
the quantum dynamics. This field theoretic model leads to the view of the phase transition as a condensation 
that is comparable to the formation of fog and rain drops from water vapor, and that might serve to model both 
the gamma and beta phase transitions. According to such a model, the production of activity with long-range 
correlation in the brain takes place through the mechanism of spontaneous breakdown of symmetry (SBS), which 
has for decades been shown to describe long-range correlation in condensed matter physics. The adoption of such 
a field theoretic approach enables modelling of the whole cerebral hemisphere and its hierarchy of components 
down to the atomic level as a fully integrated macroscopic quantum system, namely as a macroscopic system which 
is a quantum system not in the trivial sense that it is made, like all existing matter, by quantum components 
such as atoms and molecules, but in the sense that some of its macroscopic properties can best be described 
with recourse to quantum dynamics (see [21] and references therein). Also, according to Freeman and Vitielo, 
many-body quantum field theory appears to be the only existing theoretical tool capable to explain the dynamic 
origin of long-range correlations, their rapid and efficient formation and dissolution, their interim stability in 
ground states, the multiplicity of coexisting and possibly non-interfering ground states, their degree of ordering, 
and their rich textures relating to sensory and motor facets of behaviors. It is historical fact that many-body 
quantum field theory has been devised and constructed in past decades exactly to understand features like ordered 
pattern formation and phase transitions in condensed matter physics that could not be understood in classical 
physics, similar to those in the brain. 

^Our entropic action-amplitude formalism represents a kind of a generalization of the Haken-Kelso-Bunz 
(HKB) model of self-organization in the individual's motor system [23. .63^. including: multi-stability, phase 
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On the other hand, it is well-known that humans possess more degrees of freedom than 
are needed to perform any defined motor task, but are required to co-ordinate them in order 
to reliably accomplish high-level goals, while faced with intense motor variability. In an 
attempt to explain how this takes place, Todorov and Jordan have formulated an alternative 
theory of human motor co-ordination based on the concept of stochastic optimal feedback 
control |80| . They were able to conciliate the requirement of goal achievement (e.g., grasping 
an object) with that of motor variability (biomechanical degrees of freedom). Moreover, 
their theory accommodates the idea that the human motor control mechanism uses internal 
'functional synergies' to regulate task-irrelevant (redundant) movement. 

Also, a developing field in coordination dynamics involves the theory of social coordina- 
tion, which attempts to relate the DC to normal human development of complex social cues 
following certain patterns of interaction. This work is aimed at understanding how human 
social interaction is mediated by meta-stability of neural networks. fMRI and EEG are par- 
ticularly useful in mapping thalamocortical response to social cues in experimental studies. 
In particular, a new theory called the Phi complex has been developed by S. Kelso and 
collaborators, to provide experimental results for the theory of social coordination dynamics 
(see the recent nonlinear dynamics paper discussing social coordination and EEG dynamics 
|81p. According to this theory, a pair of phi rhythms, likely generated in the mirror neuron 
system, is the hallmark of human social coordination. Using a dual-EEG recording system, 
the authors monitored the interactions of eight pairs of subjects as they moved their fingers 
with and without a view of the other individual in the pair. 

Finally, the chaotic psycho-physical phase transitions embedded in CD may give a formal 
description for a phenomenon called crowd turbulence by D. Helbing, depicting crowd disas- 
ters caused by the panic stampede that can occur at high pedestrian densities and which is 
a serious concern during mass events like soccer championship games or annual pilgrimage 
in Makkah (see [35l [36l EH [60] ) . 

transitions and hysteresis effects, presenting a contrary view to the purely feedback driven systems. HKB uses 
the concepts of synergetics (order parameters, control parameters, instability, etc) and the mathematical tools of 
nonlinearly coupled (nonlinear) dynamical systems to account for self-organized behavior both at the cooperative, 
coordinative level and at the level of the individual coordinating elements. The HKB model stands as a building 
block upon which numerous extensions and elaborations have been constructed. In particular, it has been possible 
to derive it from a realistic model of the cortical sheet in which neural areas undergo a reorganization that is 
mediated by intra- and inter-cortical connections. Also, the HKB model describes phase transitions ('switches') in 
coordinated human movement as follows: (i) when the agent begins in the anti-phase mode and speed of movement 
is increased, a spontaneous switch to symmetrical, in-phase movement occurs; (ii) this transition happens swiftly 
at a certain critical frequency; (iii) after the switch has occurred and the movement rate is now decreased the 
subject remains in the symmetrical mode, i.e. she does not switch back; and (iv) no such transitions occur if the 
subject begins with symmetrical, in-phase movements. The HKB dynamics of the order parameter relative phase 
as is given by a nonlinear first-order ODE: 

(j> = {a + 2f3r^) sin — /3r^ sin 20, 

where is the phase relation (that characterizes the observed patterns of behavior, changes abruptly at the 
transition and is only weakly dependent on parameters outside the phase transition), r is the oscillator amplitude, 
while a, l3 are coupling parameters (from which the critical frequency where the phase transition occurs can be 
calculated) . 
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3 Formal crowd dynamics 

In this section we formally develop a three-step crowd psycho-physical dynamics, concep- 
tualized by transition maps in agreement with Haken's synergetics [Ml Hi]- We 
first develop a macro-level individual psycho-physical dynamics XT). Then we generalize 
I'D into an 'orchestrated' behavioral-compositional crowd dynamics CP, using a quantum- 
like micro-level formalism with individual agents representing 'crowd quanta'. Finally we 
develop a meso-level aggregate statistical-field dynamics AD, such that composition of the 
aggregates AV makes-up the crowd. 

3.1 Individual behavioral dynamics {I'D) 

IT) transition map ^ is developed using the following action-amplitude formalism (see 



1. Macroscopically, as a smooth Riemannian n— manifold Mid (see Appendix) with steady 
force-fields and behavioral paths, modelled by a real-valued classical action functional 

S\T> [^] , of the form 



(where macroscopic paths, fields and geometries are commonly denoted by an abstract 
field symbol with the potential-energy based Lagrangian L given by 



where C is Lagrangian density, the integral is taken over all n local coordinates = 
x^{t) of the ID, and d^j^^ are time and space partial derivatives of the variables 
over coordinates. The standard least action principle 



gives, in the form of the Euler-Lagrangian equations, a shortest path, an extreme force- 
field, with a geometry of minimal curvature and topology without holes. We will see 
below that high Riemannian curvature generates chaotic behavior, while holes in the 
manifold produce topologically induced phase transitions. 
2. Microscopically, as a collection of wildly fiuctuating and jumping paths (histories), 
force-fields and geometries/topologies, modelled by a complex-valued adaptive path 
integral, formulated by defining a multi-phase and multi-path (multi-field and multi- 
geometry) transition amplitude from the entropy-growing state of Mental Preparation 
to the entropy-conserving state of Physical Action, 



ElISI]): 






(6) 



where the functional ID-measure I?[i(;<I>] is defined as a weighted product 

N 

V[w^] = lim I I Wsd^l, {i — 1, n — con + dis), 



(7) 
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representing an oo— dimensional neural network |52| . with weights Ws updating by the 
general rule 

new value{t + I) ~ old value(t) + innovation(t). 

More precisely, the weights Wg = Wg (t) in ^ are updated according to one of the two 
standard neural learning schemes, in which the micro-time level is traversed in discrete 
steps, i.e., if t = tQ,ti, ...,ts then t + 1 = ti,t2, ts+iEl 

(a) A self-organized, unsupervised (e.g., Hebbian-like |34j ) learning rule: 

Ws{t + l) = wg{t) + ^{wt{t)-w:(t)), (8) 

where a = cr(t), rj = ri{t) denote signal and noise, respectively, while superscripts 
d and a denote desired and achieved micro-states, respectively; or 

(b) A certain form of a supervised gradient descent learning: 

Ws{t + l) = Ws{t)-r^'sJJ{t), (9) 

where is a small constant, called the step size, or the learning rate, and VJ(n) 
denotes the gradient of the 'performance hyper-surface' at the t—ih iteration. 

(Note that we could also use a reward-based, reinforcement learning rule 79J, in which 
system learns its optimal policy: innovation{t) = \reward{t) — penalty {t)\. ) 

In this way, we effectively derive a unique and globally smooth, causal and entropic phase- 
transition map performed at a macroscopic (global) time-level from some initial time 
tini to the final time t/i„. Thus, we have obtained macro-objects in the ID: a single path de- 
scribed by Newtonian-like equation of motion, a single force-field described by Maxwellian- 
like field equations, and a single obstacle-free Riemannian geometry (with global topology 
without holes). 

In particular, on the macro-level, we have the ID-paths, that is biodynamical trajectories 
generated by the Hamilton action principle 

(55idH = 0, 

with the Newtonian action Siu [x] given by (Einstein's summation convention over repeated 
indices is always assumed) 

/•*/'" 1 . . 
Siu[x] = y + -g^J x'x^] dt, (10) 

where ip — if{t, x*) denotes the mental LSF-potential field, while the second term, 

T 2 ' 

*The traditional neural networks approaches are known for their classes of functions they can represent. Here 
we are talking about functions in an extensional rather than merely intensional sense; that is, function can be 
read as input/output behavior ^[5l 1181 133] . This limitation has been attributed to their low-dimensionality (the 
largest neural networks are limited to the order of 10^ dimensions [59]). The proposed path integral approach 
represents a new family of function-representation methods, which potentially offers a basis for a fundamentally 
more expansive solution. 
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represents the physical (biodynamic) kinetic energy generated by the Riemannian inertial 
metric tensor gij of the configuration biodynamic manifold Mid (see Figure [1]). The corre- 
sponding Euler-Lagrangian equations give the Newtonian equations of human movement 

jT^.-T,.=F,, (11) 

where subscripts denote the partial derivatives and we have defined the covariant muscular 
forces Fi — Fi(t,x'^,x^) as negative gradients of the mental potential (p{x^), 

F^-^,.. (12) 

Equation (jlip can be put into the standard Lagrangian form as 

^Li. =L,., with L = T-^ix'), (13) 

or (using the Legendre transform) into the forced, dissipative Hamiltonian form [42l |45] 

^dp^H + dp^ R, Pi=F,- 9^. H + d^.R, (14) 

where Pi are the generalized momenta (canonically-conjugate to the coordinates x'), H = 
H{p, x) is the Hamiltonian (total energy function) and R = R{p, x) is the general dissipative 
function. 

The human motor system possesses many independently controllable components that 
often allow for more than a single movement pattern to be performed in order to achieve 
a goal. Hence, the motor system is endowed with a high level of adaptability to different 
tasks and also environmental contexts |40j. The multiple SE(3)-dynamics applied to human 
musculo-skeletal system gives the fundamental law of biodynamics, which is the covariant 
force law. 

Force co- vector field — Mass distribution x Acceleration vector- field, (15) 

which is formally written: 

Fi = gija\ {ij = 1, = dim(Af)) 

where Fi are the covariant force/torque components, gij is the inertial metric tensor of 
the configuration Riemannian manifold M = Yii SE{3y {gij defines the mass-distribution 
of the human body), while are the contravariant components of the linear and angular 
acceleration vector-field. (This fundamental biodynamic law states that contrary to common 
perception, acceleration and force are not quantities of the same nature: while acceleration 
is a non-inertial vector-field, force is an inertial co-vector-field. This apparently insignificant 
difference becomes crucial in injury prediction/prevention, especially in its derivative form in 
which the 'massless jerk' (= a) is relatively benign, while the 'massive jolt' (= F) is deadly.) 
Both Lagrangian and (topologically equivalent) Hamiltonian development of the covariant 
force law is fully elaborated in [45l |46l [471 [50] . This is consistent with the postulation that 
human action is guided primarily by natural law [64j . 
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S£(2) S£(2) 
S£(3} SE(3) 

Figure 1: Riemannian configuration manifold Mid of human biodynamics is defined as a topo- 
logical product M = ]^ • SE{3y of constrained Euclidean 5£'(3)-groups of rigid body motion in 
3D Euclidean space (see [471 150)). acting in all major (synovial) human joints. The manifold M 
is a dynamical structure activated/controlled by potential covariant forces (fT2]) produced by a 
synergetic action of about 640 skeletal muscles [35] . 




On the micro-ID level, instead of each single trajectory defined by the Newtonian equa- 
tion of motion (jlip . we have an ensemble of fluctuating and crossing paths on the config- 
uration manifold M with weighted probabilities (of the unit total sum). This ensemble of 
micro-paths is defined by the simplest instance of our adaptive path integral ([S]), similar to 
the Feynman's original sum over histories, 

(Physical Action | Mental Preparation) m ^ V[wx] e'-^'^'l , (16) 

Jib 

where is the functional ID-measure on the space of all weighted paths, and the expo- 

nential depends on the action 5*10 [a^] given by (jlOp . 

3.2 Crowd behavioral— compositional dynamics (CV) 

In this subsection we develop a generic crowd CP, as a unique and globally smooth, causal 
and entropic phase-transition map ((S]), in which agents (or, crowd's individual entities) can 
be both humans and robots. This crowd psycho-physical action takes place in a crowd 
smooth Riemannian 3n— manifold M. Recall from Figure [T] that each individual segment of 
a human body moves in the Euclidean 3-space R'^ according to its own constrained SE(3)- 
group. Similarly, each individual agent's trajectory, = a;*(t), i = l,...n, is governed by 
the Euclidean SE(2)-group of rigid body motions in the plane. (Recall that a Lie group 
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SE(2) = SO{2) X M is a set of all 3x3- matrices of the form: 

cos 9 sin 6 x 

— sin 9 cos 9 y 

1 

including both rigid translations (i.e., Cartesian x, coordinates) and rotation matrix 



cos 9 sin ( 
— sin 9 cos ( 



in Euclidean plane (see [47l [50]).) The crowd configuration manifold M is defined as a 
union of Euclidean SE(2)-groups for all n individual agents in the crowd, that is crowd's 
configuration 3n— manifold is defined as a set 

n n 

M = ^ S'£;(2)'= = ^ 50(2)'= X R^ (17) 
fc=i fe=i 
coordinated by x'' = {x'' ,y'' ,9''}, {ior k ^ 1,2, ...,n). 

In other words, the crowd configuration manifold M is a dynamical planar graph with individ- 
ual agents' SE(2)-groups of motion in the vertices and time-dependent inter- agent distances 
lij = [x''{ti) — x^{tj)] as edges. 

Similarly to the individual case, the crowd action functional includes mental cognitive 
potential and physical kinetic energy, formally given by (with i, j = 1, 3n): 

A[x\x^;k,tj] = ]^ I I 5{ll) x\U)x^{tj) dUdt, + \ 1 9.,x'{t)x={t)dt, (18) 



Hi Jti ^ Jt 

^ - — I T'lt - 1 — It ■ 1^ ^ 



with /2 = {x^ {U ) - x^ {tj )] ^ where IN < t, ,tj,t< OUT. 



The first term in (|18[) represents the mental potential for the interaction between any two 
agents a;* and a;* within the total crowd matrix x^^ . (Although, formally, this term contains 
cognitive velocities, it still represents 'potential energy' from the physical point of view.) It 
is defined as a double integral over a delta function of the square of interval between two 
points on the paths in their individual cognitive LSFs. Interaction occurs only when this 
LSF-distance between the two agents x^ and x^ vanishes. Note that the cognitive intentions 
of any two agents generally occur at different times and tj unless ti = tj, when cognitive 
synchronization occurs. This term effectively represents the crowd cognitive controller (see 

ED). 

The second term in represents kinetic energy of the physical interaction of agents. 
Namely, after the above cognitive synchronization is completed, the second term of physical 
kinetic energy is activated in the common CD manifold, reducing it to just one of the 
agents' individual manifolds, which is equivalent to the centcr-of-mass segment in the human 
musculo-skeletal system. Therefore, from (|18p we can derive a generic Euler-Lagrangian 
dynamics that is a composition of ()13|) . which also means that we have in place a generic 
Hamiltonian dynamics that is a amalgamate of (|14p , as well as the crowd covariant force law 
IjlSp . the governing law of crowd biodynamics: 

Crowd force co- vector field — Crowd mass distribution x Crowd acceleration vector-field, 
formally: Fi = gija-' , where gij is the inertial metric tensor of crowd manifold (M)) 
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The left-hand side of this equation defines forces acting on the crowd, while right-hand defines 
its mass distribution coupled to the crowd kinematics {CIC, described in the next subsection). 

At the slave level, the adaptive path integral, representing an oo— dimensional neural 
network, corresponding to the psycho-physical crowd action dTH]) , reads 

(Physical Action | Mental Preparation) cd = / x, y] Q^^l^^vM^t,] ^ ^20) 

where the Lebesgue-type integration is performed over all continuous paths — x^{ti) and 
— y-'itj), while summation is performed over all associated discrete Markov fluctuations 
and jumps. The symbolic differential in the path integral (|20p represents an adaptive path 
measure, defined as the weighted product 

N 

T>[w, X, y] ^ lim Y[w-jdx'dy^ , (i, j = 1, n). (21) 

s=l 

The quantum-field path integral (P0)) - (PT|) defines the microstate CP— level, an ensemble of 
fluctuating and crossing paths on the crowd Sri— manifold M . 

The crowd manifold M itself has quite a sophisticated topological structure defined by its 
macrostate Euler-Lagrangian dynamics. As a Riemannian smooth n— manifold, M gives rise 
to its fundamental n— groupoid, or n— category n„(M) (see [47l[50]). In n„(M), 0-cells are 
points in M ; 1-cells are paths in M (i.e., parameterized smooth maps / : [0, 1] M); 2~cells 
are smooth homotopies (denoted by ~) of paths relative to endpoints (i.e., parameterized 
smooth maps h : [0, 1] x [0, 1] M); 3-cells are smooth homotopies of homotopies of paths in 
M (i.e., parameterized smooth maps j ; [0, 1] x [0, 1] x [0, 1] M). Categorical composition 
is defined by pasting paths and homotopies. In this way, the following recursive homotopy 
dynamics emerges on the crowd 3n— manifold M: 
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/ 

2 — cell : xq • h • xi h : f ^ g d M, 

9 



h : [0, 1] X [0, 1] ^ M, : / ^ 5, 5 = h{f{xo)), 
h{xo, 0) = /(a;o), ^(a^o, 1) = 5(2^0), ^(0, t) = xq, h{l, t) = xi 
e.g., linear homotopy: h{xo,t) — (1 —t)f{xo) +tg{xQ); or 
homotopy between two Euler-Lagrangian (/, g) — dynamics 
with the same endpoint conditions (xo ,xi): 

4-/i'=/x^ and ^g±i = gxi with x{0) ^ xq, = xi; 



J : [0, 1] X [0, 1] X [0, 1] ^ M, j:h^t,i= j{h{f{xo))) 
j{xa,t,0) = h{f{xo)), j{xo,t, 1) i(/(xo)), 
j{xo,0,s) = f{xo), j{xo,l,s) = g{xo), 

j{0,t,s) = Xq, j{l,t,s) = Xi 

e.g., hnear composite homotopy: j{xo,t,s) = (1 — t)h{f{xo)) +ti{f{xo)); 
or, homotopy between two homotopies between above two Euler- 
Lagrangian (/,(?) — dynamics with the same endpoint conditions (a;o,a;i). 

3.3 Dissipative crowd kinematics (C/C) 

The crowd action (fT8|) with its amalgamate Lagrangian dynamics ([13]) and amalgamate 
Hamiltonian dynamics (HU, as well as the crowd force law define the macroscopic crowd 



/ 
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dynamics, CD. Suppose, for a moment, that CD is force-free and dissipation free, there- 
fore conservative. Now, the basic characteristic of the conservative Lagrangian/Hamiltonian 
systems evolving in the phase space spanned by the system coordinates and their veloci- 
ties/momenta, is that their flow (pf' (explained below) preserves the phase-space volume, 
as proposed by the Liouville theorem, which is the well-known fact in statistical mechanics. 
However, the preservation of the phase volume causes structural instability of the conser- 
vative system, i.e., the phase-space spreading effect by which small phase regions Rt will 
tend to get distorted from the initial one Ro during the conservative system evolution. This 
problem, governed by entropy growth {dtS > 0), is much more serious in higher dimensions 
than in lower dimensions, since there are so many 'directions' in which the region can locally 
spread (see [701 147 ) ). This phenomenon is related to conservative Hamiltonian chaos (see 
section |4] below) . 

However, this situation is not very frequent in case of 'organized' human crowd. Its self- 
organization mechanisms are clearly much stronger than the conservative statistical mechan- 
ics effects, which we interpret in terms of Prigogine's dissipative structures (see Appendix). 
Formally, if dissipation of energy in a system is much stronger then its inertial characteristics, 
then instead of the second-order Newton-Lagrangian dynamic equations of motion, we are 
actually dealing with the first-order driftless (non-acceleration, non-inertial) kinematic equa- 
tions of motion (see Appendix, eq. (|5D|) ). which is related to dissipative chaos ^7\. Briefly, 
the dissipative crowd flow can be depicted like this: from the set of initial conditions for 
individual agents, the crowd evolves in time towards the set of the corresponding entangled 
attractorsl^ which are mutually separated by fractal (non-integer dimension) separatrices. 

^Recall that quantum entanglement is a quantum mechanical phenomenon in which the quantum states of two 
or more objects are linked together so that one object can no longer be adequately described without full mention 
of its counterpart - even though the individual objects may be spatially separated. This interconnection leads to 
correlations between observable physical properties of remote systems. The related phenomenon of wave-function 
collapse gives an impression that measurements performed on one system instantaneously influence the other 
systems entangled with the measured system, even when far apart. 

Entanglement has many applications in quantum information theory. Mixed state entanglement can be viewed 
as a resource for quantum communication. A common measure of entanglement is the entropy of a mixed quantum 
state (see, e.g. [57] )• Since a mixed quantum state p is a probability distribution over a quantum ensemble, this 
leads naturally to the definition of the von Neumann entropy, S{p) — — Tr(plog2p), which is obviously 
similar to the classical Shannon entropy for probability distributions (pi,-- - ,Pn), defined as S{pi,--- ,pn) ~ 
— ^jPilogjPi. As in statistical mechanics, one can say that the more uncertainty (number of microstates) the 
system should possess, the larger is its entropy. Entropy gives a tool which can be used to quantify entanglement. 
If the overall system is pure, the entropy of one subsystem can be used to measure its degree of entanglement 
with the other subsystems. 

The most popular issue in a research on dissipative quantum brain modelling has been quantum entanglement 
between the brain and its environment [731 174) . where the brain-environment system has an entangled 'memory' 
state, identified with the ground (vacuum) state |0 >Ar, that cannot be factorized into two single-mode states. 
(In the Vitiello-Pessa dissipative quantum brain model [731 174[ , the evolution of the A/'-coded memory system 
was represented as a trajectory of given initial condition running over time-dependent states |0(t) >Af, each one 
minimizing the free energy functional.) Similar to this microscopic brain-environment entanglement, we propose 
a kind of macroscopic entanglement between the operating modes of the crowd psycho-physical controller and its 
biodynamics, which can be considered as a 'long-range correlation'. 

Applied externally to the dimension of the crowd 3n— manifold M, entanglement effectively reduces the number 
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In this subsection we elaborate on the dissipative crowd kinematics (CJC), which is self- 
controlled and dominates the CD if the crowd's inertial forces are much weaker then the the 
crowd's dissipation of energy, presented here in the form of nonlinear velocity controllers. 

Recall that the essential concept in dynamical systems theory is the notion of a vector- 
field (that we will denote by a boldface symbol), which assigns a tangent vector to each point 
p in the manifold in case. In particular, v is a gradient vector-field if it equals the gradient 
of some scalar function. A flow-line of a vector-field v is a path "f(t) satisfying the vector 
ODE, ^{t) = v{'y{t)), that is, v yields the velocity field of the path j{t). The set of all 
flow hnes of a vector-field v comprises its flow (pt that is (technically, see e.g., [171 [50]) 
a one-parameter Lie group of diffeomorphisms (smooth bijective functions) generated by a 
vector- field v on Af, such that 

ft°fs = ft+s, V'o = identity, which gives: j{t) = (pt{j{0)). 

Analytically, a vector-field v is defined as a set of autonomous ODEs. Its solution gives the 
flow ft, consisting of integral curves (or, flow hnes) ■^{t) of the vector-field, such that all the 
vectors from the vector-field are tangent to integral curves at different representative points 
p G M. In this way, through every representative point p d M passes both a curve from the 
flow and its tangent vector from the vector-field. Geometrically, vector-field is defined as a 
cross-section of the tangent bundle TM of the manifold M. 

In general, given an nD frame {di} = {d/dx''} on a smooth n— manifold M (that is, a 
basis of tangent vectors in a local coordinate chart = {x^, x") C M), we can define any 
vector- field v on M by its components — w*(t) as 

* dx"^ dx^ ^ ^ dx"' 

Thus, a vector-field v € X{M) (where X{M) is the set of all smooth vector-fields on M) is 
actually a differential operator that can be used to differentiate any smooth scalar function 
/ = /(^^i ...,x") on M, as a directional derivative of / in the direction of v. This is denoted 
simply v/, such that 

= V Oif = V — = v'—r + ... -fv"— — . 

ox^ ox^ ax^ 

In particular, if v = j{t) is a velocity vector-field of a space curve 7(t) = {x^ (t), x"{t)), 
defined by its components v'^ = x'^{t), directional derivative of /(x*) in the direction of v 
becomes 

which is a rate-of-change of / along the curve ^{t) at a point x^{t). 

Given two vector-fields, u = u^di,v = v^di G X{M), their Lie bracket (or, commutator) 
is another vector-field [u, v] e X{M), defined by 

[u, v] = uv — vu = u^di v-'dj — v^dj u^di, 



of active degrees of freedom in (|17p . 
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which, apphed to any smooth function / on Af, gives 

[u,v](/)=u(v(/))-v(u(/)). 

The Lie bracket measures the failure of 'mixed directional derivatives' to commute. 
Clearly, mixed partial derivatives do commute, [di,dj\ = 0, while in general it is not the 
case, [u,v] ^ 0. In addition, suppose that u generates the flow Lpt and v generates the flow 
LPs- Then, for any smooth function / on M, we have at any point p on M, 

[u,v](/)(p) - ^ - f{^t{^M% 

which means that in f{^s{^t{p)) we are starting at p, flowing along v a little bit, then along 
u a little bit, and then evaluating /, while in f{ipt{ips{p)) we are flowing first along u and 
then V. Therefore, the Lie bracket infinitesimally measures how these flows fail to commute. 

The Lie bracket satisfies the following three properties (for any three vector-fields u, v, w G 
M and two constants a,b - thus forming a Lie algebra on the crowd manifold M): 

(i) [u, v] = — [v, u]— skew-symmetry; 

(ii) [u, av + 6w] = a[u, v] -|- b[u, w]— bilinearity; and 

(iii) [u, [v, w]] + [v, [w, u]] + [w, [u, v]]— Jacobi identity. 

A new set of vector-fields on M can be generated by repeated Lie brackets of u, v, w G M. 

The Lie bracket is a standard tool in geometric nonlinear control theory (see, e.g. [171 ISU] ). 
Its action on vector-fields can be best visualized using the popular car parking example, in 
which the driver has two different vector-field transformations at his disposal. They can 
turn the steering wheel, or they can drive the car forward or backward. Here, we specify the 
state of a car by four coordinates: the {x,y) coordinates of the center of the rear axle, the 
direction 6 of the car, and the angle between the front wheels and the direction of the car. 
/ is the constant length of the car. Therefore, the 4D configuration manifold of a car is a set 
M = SO{2) X R^, coordinated by x = {x,y,6,(j)}, which is slightly more complicated than 
the individual crowd agent's 3D configuration manifold SE{2) = 5*0(2) x R, coordinated by 
X = {x, y, 9}. The driftless car kinematics can be defined as a vector ODE: 

x = u(x)ci +v(x)c2, (22) 

with two vector-fields, u, v G X{M), and two scalar control inputs, ci and C2. The infinites- 
imal car-parking transformations will be the following vector-fields 

\ 
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, , n ■ n tan^ d 

u(x) = DRIVE = cos 6*— -I- sm 6*— H ; — — 

ox ay I oO 



cos f 
sin 6 
7 tanc 



and 



v(x) 



STEER 



d_ 

d<j) 






V 1 / 



The car kinematics (p2)) therefore expands into a matrix ODE: 



/ i \ 

y 
e 

\4> I 



= DRIVE • Cl 





/ COS 61 \ 




f M 




sin 









J tan 


•Cl + 







I ) 




1 1 / 



• C2 



However, steer and drive do not commute (otherwise we could do aU your steering at 
home before driving of on a trip). Their combination is given by the Lie bracket 



[v, u] = [steer, drive] 



d 



I cos^ (h do 



WRIGGLE. 



The operation [v, u] = wriggle = [steer, drive] is the infinitesimal version of the sequence 
of transformations: steer, drive, steer back, and drive back, i.e., 

{steer, drive, steer^^, drive^^}. 

Now, WRIGGLE can get us out of some parking spaces, but not tight ones: we may not 
have enough room to wriggle out. The usual tight parking space restricts the DRIVE 
transformation, but not STEER. A truly tight parking space restricts steer as well by 
putting your front wheels against the curb. 

Fortunately, there is still another commutator available: 



[U: [V, U]] 

[drive, wriggle] 



[drive, [steer, drive]] = [[u, v],u] 

1 



I cos^ 



sm^- cosO— I = SLIDE. 

ax ay 



The operation [[u, v],u] = slide = [drive, wriggle] is a displacement at right angles to 
the car, and can get us out of any parking place. We just need to remember to steer, drive, 
steer back, drive some more, steer, drive back, steer back, and drive back: 

{steer, drive, steers ^, drive, steer, drive^^, steer^\ drive^^}. 

We have to reverse steer in the middle of the parking place. This is not intuitive, and no 
doubt is part of a common problem with parallel parking. 

Thus, from only two controls, ci and C2, we can form the vector-fields drive = u, steer 
= V, wriggle = [v, u], and SLIDE = [[u,v],u], allowing us to move anywhere in the car 
configuration manifold M = SO{2) x E^. All above computations are straightforward in 
M athematica^ if we define the following three symbolic functions: 

1. Jacobian matrix: JacMat[v_List, x_List] := Outer [D, v, x]; 

2. Lie bracket: LieBrc[uXist, vXist, x_List] := JacMat[v, x] . u - JacMat[u, x] . v; 

3. Repeated Lie bracket: Adj[u_List, v_List, xXist, k_] := 

If[k == 0, V, LieBrc[u, Adj[u, v, x, k - 1], x] ]; 



The above computations could instead be done in other available packages, such as Maple, by suitably 
translating the provided example code. 



19 



In case of the human crowd, we have a sHghtly simpler, but multiphed problem, i.e., 
superposition of n individual agents' motions. So, we can define the dissipative crowd kine- 
matics as a system of n vector ODEs: 

= u'' (x) c'l + v'' (x) 4 , where (23) 



d . d 



cos'' I 



u'''(x) = drive'' = cos'' 6'-— r + sin''^ ^— -r = sin'' 6* I, and 

v''(x) = steer'' = = I I , while c\ and c\ are crowd controls. 
Thus, the crowd kinematics ((23|) expands into the matrix ODE: 

X 

y 






1 cos'' 6* ' 










drive'' • Ci - 


f steer'' -02= sin'' e 




(•) 




(24) 




I , 











The dissipative crowd kinematics obeys the set of n-tuple integral rules of 

motion that are similar (though slightly simpler) to the above rules of the car kinematics, 
including the following derived vector-fields: 

wriggle''' = [steer'', drive'''] = [v''',u''] and slide'' = [drive''', wriggle''] = [[u'',v''],i 

Thus, controlled by the two vector controls and C2, the crowd can form the vector- 
fields: drive = u'', steer = v'', WRIGGLE = [v'',u''], and slide = [[u'', v''], u''], allowing 
it to move anywhere within its configuration manifold M given by (|17[) . Solution of the 
dissipative crowd kinematics ([23 |) -([24 |) defines the dissipative crowd fiow, . 

Now, the general CD-CK. crowd behavior can be defined as a amalgamate fiow (psycho- 
physical Lagrangian fiow, 0f , plus dissipative kinematic flow, (f)^) on the crowd manifold M 
defined by dTT]), 

Ct^4>t +4>t ■t^{M{t),g{t)), 

which is a one-parameter family of homcomorphic (topologically equivalent) Riemannian 
manifoldf0 {M,g = gtj), parameterized by a 'time' parameter t. That is, Ct can be used for 



^Proper differentiation of vector and tensor fields on a smooth Riemannian manifold (like the crowd 
3n— manifold M) is performed using the Levi-Civita covariant derivative (see, e.g., [471150) ). Formally, let M be a 
Riemannian A^— manifold with the tangent bundle TM and a local coordinate system {a;*}^! defined in an open 
set U C M. The covariant derivative operator, Vx '■ C°°{TM) — > C°°{TM), is the unique linear map such that 
for any vector-fields X, Y, Z, constant c, and scalar function / the following properties are valid: 

Vx+oY = \7x +cVy, \7x{Y + fZ)^VxY + {Xf)Z + fVxZ, VxY ^VyX ^ [X,Y], 

where [X, Y] is the Lie bracket of X and Y. In local coordinates, the metric g is defined for any orthonormal basis 
{di = d/dx^) in U C M by gij = g{di,dj) = Sij, dugij ~ 0. Then the affine Levi-Civita connection is defined on 
M by 

Va^dj = Fijdk, where fJj = — g''' (digji + djgu — digij) are the Christoffel symbols. 
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describing smooth deformations of the crowd manifold M over time. The manifold family 
{M{t), g{t)) at time t determines the manifold family {M{t + dt), g{t + dt)) at an infinitesimal 
time t + dt into the future, according to some presecribed geometric flow, like the celebrated 
Ricci flow [221 [301 1321 131] (that was an instrument for a proof of a 100-year old Poincare 
conjecture) , 

^tg^J{t) = -2i?y(i), (25) 

where Rij is the Ricci curvature tensor (see Appendix) of the crowd manifold M and dtg(t) 
is defined as 

dtg{t) = —g{t) := hm . 26 

dt dt^o dt 

3.4 Aggregate behavioral— compositional dynamics {AV) 

To formally develop the meso- level aggregate behavioral-compositional dynamics (-42?), we 
start with the crowd path integral (|20p . which can be redefined if we Wick-rotate the time 
variable t to imaginary values, t ^ t — it, thereby transforming the Lorentzian path integral 
in real time into the Euclidean path integral in imaginary time. Furthermore, if we rectify the 
time axis back to the real line, we get the adaptive SFT-partition function as our proposed 
^P-model: 

(Physical Action | Mental Preparation) ad = / Vlw, x, y] e-^[^'f'*"*jl. (27) 

JcD 



Now, using the covariant derivative operator Vjc we can define the Riemann curvature (3, 1)— tensor JRm by 

D\m{X,Y)Z — Vx^yZ ~ Vy^ xZ — V[x,Y]Z, 

which measures the curvature of the manifold by expressing how noncommutative covariant differentiation is. 
The (3, 1)— components liij}^ of SHm are defined in (7 C M by 

d\m{dt,dj)dk = Rljkdi, or R\jk = diV]^ - djV\k +V'Jl:V\^ -V'S.T^rn- 
Also, the Riemann (4, 0)— tensor Rijki ~ gimRTjk is defined as the g— based inner product on M, 

Rijki = (D^m {di, dj)dk,di) . 
The first and second Bianchi identities for the Riemann (4,0)— tensor Rijki hold, 

Rijki + Rjkil + Rkijl ~ 0, ^iRjklm + ^ j Rkilm + ^ kRijlm = 0, 

while the twice contracted second Bianchi identity reads: 2\/jRij = ViR. 
The (0, 2) Ricci tensor Die is the trace of the Riemann (3, 1)— tensor [ftm, 

%c{Y,Z)+tr{X ^^Rm{X,Y)Z), so that y\c{X,Y) = g{mm{^^, X)^^,Y), 

Its components Rjk ~ SHc {dj, dk) are given in [/ C M by the contraction 

D — R* r>r R 1^ V' — fl M V' _ P"* 

J^jk — ^ijk 5 J^jk — t^z-L jk ^fc J- ji ^ J- jni^ jk ^ mk^ ji • 

Finally, the scalar curvature R is the trace of the Ricci tensor JRc, given in [/ C M by; R — g^-' Rij. 
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The adaptive ^2?-transition amplitude (Physical Action | Mental Preparation) ad as de- 
fined by the SFT-partition function (j27p is a general model of aMarkov stochastic process. 
Recall that Markov process is a random process characterized by a lack of memory, i.e., 
the statistical properties of the immediate future are uniquely determined by the present, 
regardless of the past (see, e.g. [551137]). The A^— dimensional Markov process can be defined 
by the Ito stochastic differential equation, 

dxi(t) = A,[x'{t),t]dt + B,j[x'{t),t]dW^t), (28) 
x'iO) = x,o, (z,j = l,...,iV) (29) 

or corresponding Ito stochastic integral equation 

x\t)^x'{0)+ f dsA,[x\s),s]+ f dW^is)B,j[x'{s),s], (30) 



in which x^{t) is the variable of interest, the vector Ai[x{t),t] denotes deterministic drift, the 
matrix Bij[x{t),t] represents continuous stochastic diffusion fluctuations, and (t) is an 
N— variable Wiener process (i.e., generalized Brownian motion [22^) and 

dW^{t) = W^{t + dt) - W^{t). 

The two Ito equations (|29p - (|30p are equivalent to the general Chapman-Kolmogorov 
probability equation (see equation pip below). There are three well-known special cases of 
the Chapman-Kolmogorov equation (see [55]): 

1. When both Bij [x{t),t] and W{t) are zero, i.e., in the case of pure deterministic motion, 
it reduces to the Liouville equation 

dtP{x', t'\x", t") - - E ^ {Ar[x{t),t] Fix', t'\x", t")} . 



2. When only W{t) is zero, it reduces to the Fokker-Planck equation 

_d_ 

dx 



dtP{x',t'\x",t") = -Y^l-{AMt),t]P{x\t'\x",t")} 



+ lT.d^iB,^^-(')'t]P{x\t'\x",t")}. 



3. When both Ai[x{t),t] and Bij[x{t),t] are zero, i.e., the state-space consists of integers 
only, it reduces to the Master equation of discontinuous jumps 

dtP[x',t'\x",t")= I dxW{x'\x",t)P{x',t'\x",t")- I dxW{x"\x',t)P{x',t'\x",t"). 



The Markov assumption can now be formulated in terms of the conditional probabilities 
P{x^ ,ti): if the times ti increase from right to left, the conditional probability is determined 
entirely by the knowledge of the most recent condition. Markov process is generated by a set 
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of conditional probabilities whose probability-density P — P{x' ,t'\x" ,t") evolution obeys 
the general Chapman-Kolmogorov integro- differential equation 



dtP 



Y^ — {AMt\t]P} + 



{i?.,[x(t),t] P} 



dx^dxi 



i ij 



+ 



I 



dx {W{x'\x",t) P - W(x"\x', t) P] 



including deterministic drift, diffusion fluctuations and discontinuous jumps (given respec- 
tively in the first, second and third terms on the r.h.s.). This general Chapman-Kolmogorov 
integro-differential equation PT|) . with its conditional probability density evolution, P = 
P{x' ,t'\x" ,t"), is represented by our SFT-partition function ([?7|) . 

Furthermore, discretization of the adaptive SFT-partition function (|27p gives the stan- 
dard partition function (see Appendix) 



where is the motion energy eigenvalue (reflecting each possible motivational energetic 
state), T is the temperature-like environmental control parameter, and the sum runs over 
all ID energy eigenstates (labelled by the index j). From (I3ip . we can calculate the transition 
entropy, as S = kB in Z (see the next section). 

4 Entropy, chaos and phase transitions in the crowd 
manifold 

Recall that nonequilibrium phase transitions [HI [2Z1 HH] are phenomena which bring 

about qualitative physical changes at the macroscopic level in presence of the same mi- 
croscopic forces acting among the constituents of a system. In this section we extend the 
CD formalism to incorporate both algorithmic and geometrical entropy as well as dynam- 
ical chaos [58l |48l [56j between the entropy-growing phase of Mental Preparation and the 
entropy-conserving phase of Physical Action, together with the associated topological phase 
transitions. 

4.1 Algorithmic entropy 

The Boltzmann and Shannon (hence also Gibbs entropy, which is Shannon entropy scaled by 
fcln2, where k is the Bolzmann constant) entropy definitions involve the notion of ensembles. 
Membership of microscopic states in ensembles defines the probability density function that 
underpins the entropy function; the result is that the entropy of a definite and completely 
known microscopic state is precisely zero. Bolzmann entropy defines the probabilistic model 
of the system by effectively discarding part of the information about the system, while the 
Shannon entropy is concerned with measuring the ignorance of the observer - the amount of 
missing information - about the system. 




(31) 
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Zurek proposed a new physical entropy measure that can be apphed to individual micro- 
scopic system states and does not use the ensemble structure. This is based on the notion 
of a fixed individually random object provided by Algorithmic Information Theory and Kol- 
mogorov Complexity: put simply, the randomness K{x) of a binary string x is the length in 
terms of number of bits of the smallest program p on a universal computer that can produce 

X. 

While this is the basic idea, there are some important technical details involved with 

this definition. The randomness definition uses the prefix complexity K{.) rather than the 
older Kolmogorov complexity measure C(.): the prefix complexity K{x\y) of x given y is 
the Kolmogorov complexity C(j)^{x\y) = min{p | x = (f>u{{y,p))} (with the convention that 
C(i,^{x\y) = oc if there is no such p) that is taken with respect to a reference universal 
partial recursive function that is a universal prefix function. Then the prefix complexity 
K{x) of X is just K{x\e) where e is the empty string. A partial recursive prefix function 
(f) : M ^ N is a partial recursive function such that if 4'{p) < oo and 4>{q) < oo then p is not a 
proper prefix of q: that is, we restrict the complexity definition to a set of strings (which are 
descriptions of efii'ective procedures) such that none is a proper prefix of any other. In this 
way, all effective procedure descriptions are self- delimiting: the total length of the description 
is given within the description itself. A universal prefix function (p^is a prefix function such 
that Vn € N ^„((?/, {n,p))) = (t>n{{v^p))^ where </>„ is numbered n according to some Godel 
numbering of the partial recursive functions; that is, a universal prefix function is a partial 
recursive function that simulates any partial recursive function. Here, {x, y) stands for a 
total recusive one-one mapping from NxN into N, {x\,X2, ■ ■ ■ , Xn) = {xi, {x2, ■ ■ ■ , a;„)), N is 
the set of natural numbers, and M = {0, l}*is the set of all binary strings. 

This notion of entropy circumvents the use of probability to give a concept of entropy 
that can be applied to a fully specified macroscopic state: the algorithmic randomness of the 
state is the length of the shortest possible effective description of it. To illustrate, suppose for 
the moment that the set of microscopic states is countably infinite, with each state identified 
with some natural number. It is known that the discrete version of the Gibbs entropy (and 
hence of Shannon's entropy) and the algorithmic entropy arc asymptotically consistent under 
mild assumptions. Consider a system with a countably infinite set of microscopic states X 
supporting a probability density function P(.) so that P{x) is the probability that the system 
is in microscopic state x G X. Then the Gibbs entropy is Sg{P) = — (/c In 2) Yl P{x) log P{x) 

(which is Shannon's information-theoretic entropy H(P) scaled by fcln2). Supposing that 

P(.) is recursive, then Sa(P) = (fcln2) ^ P{x)K{x) + C, where is a constant depending 

:iex 

only on the choice of the reference universal prefix function </>. Hence, as a measure of 
entropy, the function K{.) manifests the same kind of behavior as Shannon's and Gibbs 
entropy measures. 

Zurek's proposal was of a new physical entropy measure that includes contributions from 
both the randomness of a state and ignorance about it. Assume now that we have determined 
the macroscopic parameters of the system, and encode this as a string - which can always be 
converted into an equivalent binary string, which is just a natural number under a standard 
encoding. It is standard to denote the binary string and its corresponding natural number 
interchangeably; here let x be the encoded macroscopic parameters. Zurek's definition of 
algorithmic entropy of the macroscopic state is then K{x) + H^, where = SB{x)/{k In 2), 
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where Sb^x) is the Bolzmann entropy of the system constrained by x and k is Bolzmann's 
constant; the physical version of the algorithmic entropy is therefore defined as SAix) = 
(fcln2)(i^(a;) + Hx). Here -ff^; represents the level of ignorance about the microscopic state, 
given the parameter set x; it can decrease towards zero as knowledge about the state of the 
system increases, at which point the algorithmic entropy reduces to the Bolzmann entropy. 

4.2 Ricci flow and Perelman entropy— action on the crowd manifold 

Recall that the inertial metric crowd flow, Ct : i i— > {M{t), g{t)) on the crowd 3n— manifold 
p7|) is a one-parameter family of homeomorphic Riemannian manifolds {M,g), evolving by 
the Ricci flow (P5)) - (PS)) . 

Now, given a smooth scalar function m : A/ — > R on the Riemannian crowd 3ri— manifold 
M, its Laplacian operator A is locally defined as 

where Vi is the covariant derivative (or, Levi-Civita connection, see Appendix). We say that 
a smooth function u : M x [0, T) — > M, where T 6 (0, cx)], is a solution to the heat equation 
(see Appendix, eq. ((5^ ) on M if 

dtu = Au. (32) 

One of the most important properties satisfied by the heat equation is the maximum principle, 
which says that for any smooth solution to the heat equation, whatever point-wise bounds 
hold at i = also hold for t > [12 . This property exhibits the smoothing behavior of the 
heat diffusion ([32| on M. 

Closely related to the heat diffusion ([32| is the (the Fields medal winning) Perelman 
entropy-action functional, which is on a 3n— manifold M with a Riemannian metric gij and 
a (temperature-like) scalar function / given by [71| 

£ = f {R+\Wf\^)e-fdfi (33) 

where R is the scalar Riemann curvature on M, while dfi is the volume 3n— form on M, 
defined as 

d^i = ^Jdet{g^dx^ Adx'^ A ... Adx^"-. (34) 

During the Ricci flow on the crowd manifold P7)l . that is, during the inertial 

metric crowd flow, Ct : 1 1-^ {AI{t), g{t)), the Perelman entropy functional ([55]) evolves as 

dt£^2 j 1% + V,V,/pe-/d//. (35) 

Now, the crowd breathers are solitonic crowd behaviors, which could be given by localized 
periodic solutions of some nonlinear soliton PDEs, including the exactly solvable sine-Gordon 
equation and the focusing nonlinear Schrodinger equation. In particular, the time-dependent 
crowd inertial metric gij{t), evolving by the Ricci flow g{t) given by (|25p - (|26p on the crowd 
3n— manifold M is the Ricci crowd breather, if for some ti < t2 and a > the metrics agij {ti) 
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and gij{t2) differ only by a difTeomorphism; the cases a = < l,a > \ correspond to 
steady, shrinking and expanding crowd breathers, respectively. Trivial crowd breathers, for 
which the metrics gij(ti) and gij{t2) on M differ only by diffeomorphism and scaling for each 
pair of ti and t2, are the crowd Ricci solitons. Thus, if we consider the Ricci flow (P5|) - (P5)) 
as a biodynamical system on the space of Riemannian metrics modulo diffeomorphism and 
scaling, then crowd breathers and solitons correspond to periodic orbits and fixed points 
respectively. At each time the Ricci soliton metric satisfies on M an equation of the form 

m 

Rij + cgij + Vih.j + Vjhi = 0, 

where c is a number and bi is a 1-form; in particular, when hi = ^V^a for some function a 
on Af, we get a gradient Ricci soliton. 

Define \{gij) — inf £{gij, /), where infimum is taken over all smooth /, satisfying 



/ e-^(i/i = l. (36) 
Jm 



^{gij) is the lowest eigenvalue of the operator — 4A + R. Then the entropy evolution formula 
l[55|) implies that X{gij{t)) is non-decreasing in i, and moreover, if X{ti) = A(t2), then for 
t G [^1,^2] we have Rij + ViWjf = for / which minimizes £ on M [7T]. Therefore, a steady 
breather on M is necessarily a steady soliton. 

If we define the conjugate heat operator on M as 

□ * = -d/dt-A + R 

then we have the conjugate heat equation: 0*u = 0. 

The entropy functional (|33p is nondecreasing under the coupled Ricci-diffusion flow on 
M [54] 

R iVuP 

dtgij = -2i?ij, dtu = -Au + — , (37) 

where the second equation ensures u^dfj, = 1, to be preserved by the Ricci flow g{t) 

on M. If we define u — e^i, then ([37]) is equivalent to /—evolution equation on M (the 
nonlinear backward heat equation), 

dtf^-Af + \Wf\^~R, 

which instead preserves (|36|) . The coupled Ricci-diffusion flow (|37l) is the most general 
biodynamic model of the crowd reaction-diffusion processes on M. In a recent study [1] this 
general model has been implemented for modelling a generic perception-action cycle with 
applications to robot navigation in the form of a dynamical grid. 

Perelman's functional £ is analogous to negative thermodynamic entropy [71j . Recall 
(see Appendix) that thermodynamic partition function for a generic canonical ensemble at 
temperature f3~^ is given by 

Z = J e-^^(hj{E), (38) 
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where w{E) is a 'density measure', which does not depend on /3. From it, the average energy 
is given by (E) = —dplnZ, the entropy is S = P {E) + InZ, and the fluctuation is 
CT = {{E-{E)f)^dp.\nZ. 

If we now fix a closed 3n— manifold M with a probability measure m and a metric gij (r) 
that depends on the temperature r, then according to equation 

the partition function ([38|l is given by 

lnZ = y" (-/ + ^)dm. (39) 

From dnH) we get (see [71]) 

(£;) = -r^ [ {R+ I V/p - :J^) dm, ^ = - / (^(i? + IV/P) + / - n) dm, 

(t = 2t^ I + VjVo/ — T^g,;, p dm, where dm — udV, u ^ {4ttt)^^ c^^ . 

Jm 2r 

From the above formulas, we see that the fluctuation a is nonnegative; it vanishes only 
on a gradient shrinking soliton. (E) is nonnegative as well, whenever the flow exists for all 
sufficiently small r > 0. Furthermore, if the heat function u: (a) tends to a (5— function as 
r — > 0, or (b) is a limit of a sequence of partial heat functions Ui, such that each Ui tends to 
a (5— function as t ^ > 0, and ^ 0, then the entropy S is also nonnegative. In case (a), 
all the quantities (E) , S, a tend to zero as r — > 0, while in case (b), which may be interesting 
if gijir) becomes singular at r = 0, the entropy S may tend to a positive limit. 

4.3 Chaotic inter-phase in crowd dynamics induced by its Rieman- 
nian geometry change 

Recall that CD transition map ((5]) is defined by the chaotic crowd phase-transition amplitude 



/ dtS=0 

( PHYS. ACTION 



CHAOS 



dtS>0 \ f 

MENTAL PREP. \ := J V[x]e'^^'=\ 



where we expect the inter-phase chaotic behavior (see [ST]). To show that this chaotic inter- 
phase is caused by the change in Riemannian geometry of the crowd 3n— manifold M , we 
will first simplify the CT> action functional as 

A[x] ^IJ^"" [9^3 i^i' - V{^^ i)] dt, (40) 
with the associated standard Hamiltonian, corresponding to the amalgamate version of ([1 

N 

2 



^ 1 
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where pi are the SE(2)-niomenta, canonically conjugate to the individual agents' SE(2)- 
coordinates x^, {i = l,...,3n). Biodynamics of systems with action PO)) and Hamiltonian 
(I4ip are given by the set of geodesic equations [471 150] 

^ + r;..^;^^-o^ (42) 

as ds ds 

where r*^, are the Christoffel symbols of the affine Levi-Civita connection of the Riemannian 
CD manifold M (see Appendix). In this geometrical framework, the instability of the tra- 
jectories is the instability of the geodesies, and it is completely determined by the curvature 
properties of the CD manifold M according to the Jacobi equation of geodesic deviation 
071 [5D] 

whose solution J, usually called Jacobi variation field, locally measures the distance between 
nearby geodesies; D/ds stands for the covariant derivative along a geodesic and jj,^ are 
the components of the Riemann curvature tensor of the CD manifold M . 

The relevant part of the Jacobi equation (|43p is given by the tangent dynamics equation 

mm 

J^ + ffofeo^/'^O, (z,fc=l,...,3n), (44) 
where the only non- vanishing components of the curvature tensor of the CD manifold M are 

fi'ofco = d^V/dx'dx^. (45) 

The tangent dynamics equation (|44|) can be used to define Lyapunov exponents in dynam- 
ical systems given by the Riemannian action (|40p and Hamiltonian (I4ip . using the formula 

m 

Ai = hm l/2nog{Ml,[jf{t) + jf{t)yMl,[jf{Q) + jf (0)]). (46) 

t — ^oo 

Lyapunov exponents measure the strength of dynamical chaos in the crowd psycho-physical 
behavior. The sum of positive Lyapunov exponents defines the Kolmogorov-Sinai entropy 
(see Appendix). 



4.4 Crowd nonequilibrium phase transitions induced by manifold 
topology change 

Now, to relate these results to topological phase transitions within the CT> manifold M given 
by (fTT]) , recall that any two high-dimensional manifolds My and M^r have the same topology 
if they can be continuously and diffcrentiably deformed into one another, that is if they are 
diffeomorphic. Thus by topology change the 'loss of diffeomorphicity' is meant [76]. In this 
respect, the so-called topological theorem |20j says that non-analyticity is the 'shadow' of 
a more fundamental phenomenon occurring in the system's configuration manifold (in our 
case the CT> manifold): a topology change within the family of equipotential hypersurfaces 

My = {(a;\...,a;3") e M^"! V ix\ . . . , x^"") = v}, 
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where V and are the microscopic interaction potential and coordinates respectively. This 
topological approach to PTs stems from the numerical study of the dynamical counterpart 
of phase transitions, and precisely from the observation of discontinuous or cuspy patterns 
displayed by the largest Lyapunov exponent Ai at the transition energy [13] . Lyapunov 
exponents cannot be measured in laboratory experiments, at variance with thermodynamic 
observables, thus, being genuine dynamical observables they are only be estimated in numer- 
ical simulations of the microscopic dynamics. If there are critical points of V in configuration 
space, that is points Xc = [xi, . . . ,X3„] such that VV"(a;)|^^^ = 0, according to the Morse 
Lemma [35], in the neighborhood of any critical point Xc there always exists a coordinate 
system x{t) = [x^it), ...,x^"(t)] for which [T3j 

V{x) = V{x,) -xl xl+ + • • • + xl^, (47) 

where k is the index of the critical point, i.e., the number of negative eigenvalues of the 
Hessian of the potential energy V . In the neighborhood of a critical point of the CP- 
manifold M, equation (jlTj) yields the simphfied form of (|45l) . d'^V/dx^dx^ — ±(5^, giving 
J unstable directions that contribute to the exponential growth of the norm of the tangent 
vector J. 

This means that the strength of dynamical chaos within the CP-manifold M, measured 
by the largest Lyapunov exponent Ai given by l|46)). is affected by the existence of critical 
points Xc of the potential energy V{x). However, as V{x) is bounded below, it is a good 
Morse function, with no vanishing eigenvalues of its Hessian matrix. According to Morse 
theory [35] , the existence of critical points of V is associated with topology changes of the 
hypersurfaces {M^}„gR. The topology change of the {M^}t,gR at some Vc is a necessary 
condition for a phase transition to take place at the corresponding energy value [20j. The 
topology changes implied here are those described within the framework of Morse theory 
through 'attachment of handles' [38] to the CP-manifold M . 

In our path-integral language this means that suitable topology changes of equipotential 
submanifolds of the CP-manifold M can entail thermodynamic-like phase transitions [211 
[25l [26] , according to the general formula: 

(phase out I phase in) :— I P[w$] e'"^'*' . 

J top— ch 

The statistical behavior of the crowd biodynamics system with the action functional ([40]) and 
the Hamiltonian (I4ip is encompassed, in the canonical ensemble, by its partition function, 
given by the Hamiltonian path integral [50] 

^3n- / P[p]P[,x] exp{i / [p.x' ~H{p,x)]dT], (48) 

J top — ch J i 

where we have used the shorthand notation 

/ PWPH^/n^^^^- 

Jtop-ch J V ^'^ 
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The path integral ([48]l can be calculated as the partition function [19], 

3n , - ^ 3ti 



P OIL / \ ^ r 



where the last term is written using the so-cahed co-area formula [17] , and v labels the 
equipotential hypersurfaces of the CV manifold M, 

A'U = {ix\ x^") e M.^"\Vix\ . . . = v}. 

Equation (|49p shows that the relevant statistical information is contained in the canonical 
configurational partition function 



0V(x) 



Note that is decomposed, in the last term of (|49|) , into an infinite summation of geometric 
integrals, 

■JMv 

defined on the {Af„}„gR. Once the microscopic interaction potential V{x) is given, the 
configuration space of the system is automatically foliated into the family {Mt,}^gR of these 
equipotential hypersurfaces. Now, from standard statistical mechanical arguments we know 
that, at any given value of the inverse temperature (3, the larger the number 3n, the closer 
to My = are the microstates that significantly contribute to the averages, computed 
through Z^niP), of thermodynamic observables. The hypersurface is the one associated 
with 

-f3V(x) 



the average potential energy computed at a given /?. Thus, at any f3, if 3n is very large the 
effective support of the canonical measure shrinks very close to a single = . Hence, 
the basic origin of a phase transition lies in a suitable topology change of the {M^}, occurring 
at some Vc |19j . This topology change induces the singular behavior of the thermodynamic 
observables at a phase transition. It is conjectured that the counterpart of a phase transi- 
tion is a breaking of diffeomorphicity among the surfaces M^, it is appropriate to choose a 
diffeomorphism invariant to probe if and how the topology of the My changes as a function 
of V. Fortunately, such a topological invariant exists, the Euler characteristic of the crowd 
manifold M, defined by gllEO] 

3n 

x{M) = J2{-l)'^bk{M), (50) 

fe=0 

where the Betti numbers bk{M) are diffeomorphism invariants {bk are the dimensions of the 
de Rham's cohomology groups 77'=(M;R); therefore the bk are integers). This homological 
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formula can be simplified by the use of the Gauss-Bonnet theorem, that relates x(M) with 
the total Gauss-Kronecker curvature Kg of the CP-manifold M given by [SHI [55] 



5 Conclusion 

Our understanding of crowd dynamics is presently limited in important ways; in particular, 
the lack of a geometrically predictive theory of crowd behavior restricts the ability for au- 
thorities to intervene appropriately, or even to recognize when such intervention is needed. 
This is not merely an idle theoretical investigation: given increasing population sizes and 
thus increasing opportunity for the formation of large congregations of people, death and 
injury due to trampling and crushing - even within crowds that have not formed under com- 
mon malicious intent - is a growing concern among police, military and emergency services. 
This paper represents a contribution towards the understanding of crowd behavior for the 
purpose of better informing decision-makers about the dangers and likely consequences of 
different intervention strategies in particular circumstances. 

In this paper, we have proposed an entropic geometrical model of crowd dynamics, with 
dissipative kinematics, that operates across macro-, micro- and meso-levels. This proposi- 
tion is motivated by the need to explain the dynamics of crowds across these levels simul- 
taneously: we contend that only by doing this can we expect to adequately characterize the 
geometrical properties of crowds with respect to regimes of behavior and the changes of state 
that mark the boundaries between such regimes. 

In pursuing this idea, we have set aside traditional assumptions with respect to the 
separation of mind and body. Furthermore, we have attempted to transcend the long- 
running debate between contagion and convergence theories of crowd behavior with our 
multi-layered approach: rather than representing a reduction of the whole into parts or the 
emergence of the whole from the parts, our approach is build on the supposition that the 
direction of logical implication can and does flow in both directions simultaneously. We refer 
to this third alternative, which effectively unifies the other two, as behavioral composition. 

The most natural statistical descriptor is crowd entropy, which satisfies the extended sec- 
ond thermodynamics law applicable to open systems comprised of many components. Simi- 
larities between the configuration manifolds of individual (micro-level) and crowds (macro- 
level) motivate our claim that goal-directed movement operates under entropy conservation, 
while natural crowd dynamics operates under monotonically increasing entropy functions. 
Of particular interest is what happens between these distinct topological phases: the phase 
transition is marked by chaotic movement. 

We contend that backdrop gives us a basis on which we can build a geometrically pre- 
dictive model-theory of crowd psycho-physical behavior. This contrasts with previous ap- 
proaches, which are explanatory only (explanation that is really narrative in nature). We 
propose an entropy formulation of crowd dynamics as a three step process involving indi- 
vidual and collective psycho-dynamics, and - crucially - non-equilibrium phase transitions 
whereby the forces operating at the microscopic level result in geometrical change at the 
macroscopic level. Here we have incorporated both geometrical and algorithmic notions of 




where d^ is given by ([M]) . 
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entropy as well as chaos in studying the topological phase transition between the entropy 
conservation of physical action and the entropy increase of mental preparation. 



6 Appendix 

6.1 Extended second law of thermodynamics 

According to Boltzmann's interpretation of the second law of thermodynamics, there exists 
a function of the state variables, usually chosen to be the physical entropy S of the system 
that varies monotonically during the approach to the unique final state of thermodynamic 
equilibrium: 

dtS > (for any isolated system). (51) 

It is usually interpreted as a tendency to increased disorder, i.e., an irreversible trend to 
maximum disorder. The above interpretation of entropy and a second law is fairly obvious 
for systems of weakly interacting particles, to which the arguments developed by Boltzmann 
referred. 

However, according to Prigogine [66] , the above interpretation of entropy and a second 
law is fairly obvious only for systems of weakly interacting particles, to which the arguments 
developed by Boltzmann referred. On the other hand, for strongly interacting systems like 
the crowd, the above interpretation does not apply in a straightforward manner since, we 
know that for such systems there exists the possibility of evolving to more ordered states 
through the mechanism of phase transitions. 

Let us now turn to nonisolated systems (like a human crowd), which exchange en- 
ergy/matter with the environment. The entropy variation will now be the sum of two terms. 
One, entropy flux, di,S, is due to these exchanges; the other, entropy production, diS, is due 
to the phenomena going on within the system. Thus the entropy variation is 

a5=f + « (52) 

For an isolated system deS — 0, and (l52|) together with (l5T|) reduces to dS = diS > 0, the 
usual statement of the second law. But even if the system is nonisolated, diS will describe 
those (irreversible) processes that would still go on even in the absence of the flux term deS. 
We thus require the following extended form of the second law: 

dtS > (for any nonisolated system). (53) 

As long as diS is strictly positive, irreversible processes will go on continuously within the 
systemlfl Thus, diS > is equivalent to the condition of dissipativity as time irreversibility. 

* Among the most common irreversible processes contributing to dtS are chemical reactions, heat conduction, 
diffusion, viscous dissipation, and relaxation phenomena in electrically or magnetically polarized systems. For 
each of these phenomena two factors can be defined: an appropriate internal flux, Ji, denoting essentially its rate, 
and a driving force, Xi, related to the maintenance of the nonequilibrium constraint. A most remarkable feature 
is that diS becomes a bilinear form of Ji and Xi. The following table summarizes the fluxes and forces associated 
with some commonly observed irreversible phenomena (see [661 146) ) 



32 



If, on the other hand, diS reduces to zero, the process will be reversible and will merely join 
neighboring states of equilibrium through a slow variation of the flux term d^S. 

From a computational perspective, we have a related algorithmic entropy. Suppose we 
have a universal machine capable of simulating any effective procedure (i.e., a universal 
machine that can compute any computable function). There are several models to choose 
from, classically we would use a Universal Turing Machine but for technical reasons we 
are more interested in Lambda-type Calculi or Combinatory Logics. Let us describe the 
system of interest through some encoding as a combinatorial structure (classically this would 
be a binary string, but again I prefer for technical reasons Normal Forms with respect to 
alpha/beta/eta, weak, strong reduction, which are basically the Lambda-type Calculi and 
Combinatory Logic notions roughly akin to a "computational" step). In other words, we 
have states of our system now represented as sentences in some language. The entropy is 
simply the minimum effective procedure against our computational model that generates the 
description of the system state. This is a universal and absolute notion of compression of 
our data - the entropy is the strongest compression over all possible compression schemes, 
in effect. Now here is the 'magic': this minimum is absolute in the sense that it does not 
vary (except by a constant) with respect to our reference choice of machine. 



6.2 Thermodynamic partition function 



Recall that the partition function Z is a quantity that encodes the statistical properties of a 
system in thermodynamic equilibrium. It is a function of temperature and other parameters, 
such as the volume enclosing a gas. Other thermodynamic variables of the system, such as 
the total energy, free energy, entropy, and pressure, can be expressed in terms of the partition 
function or its derivatives. 

A canonical ensemble is a statistical ensemble representing a probability distribution 
of microscopic states of the system. Its probability distribution is characterized by the 
proportion pi of members of the ensemble which exhibit a measurable macroscopic state 
i, where the proportion of microscopic states for each macroscopic state i is given by the 
Boltzmann distribution. 



Pi 



Phenomenon 

Heat conduction 
Diffusion 
Viscous flow 
Chemical reaction 



Flux 

Heat flux, Jth 
Mass flux, Jd 
Pressure tensor, P 
Rate of reaction, u 



Force Rank 

grad{l/T) Vector 

-~[grad{fi/T) - F] Vector 

(l/T) gradv Tensor 

Aflinity of reaction Scalar 



In general, the fluxes are very complicated functions of the forces Xi . A particularly simple situation arises 
when their relation is linear, then we have the celebrated Onsager relations, 



Ji — LikXk 



(i, fc = 1, n) 



(54) 



in which Lik denote the set of phenomenological coejjicients. This is what happens near equilibrium where they 
are also symmetric, Lit = Lki- Note, however, that certain states far from equilibrium can still be characterized 
by a linear dependence of the form of (|54|) that occurs either accidentally or because of the presence of special 
types of regulatory processes. 
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where Ei is the energy of state i. It can be shown that this is the distribution which is most 

likely, if each system in the ensemble can exchange energy with a heat bath, or alternatively 
with a large number of similar systems. In other words, it is the distribution which has 
maximum entropy for a given average energy < Ei >. 

The partition function of a canonical cnsemMc is defined as a sum Z{l3) = Q-PEj ^ ■y^rl^gJ■g 
(3 = l/{kBT) is the 'inverse temperature', where T is an ordinary temperature and fcs is 
the Boltzmann's constant. However, as the position and momentum pi variables of an ith 
particle in a system can vary continuously, the set of microstates is actually uncountable. 
In this case, some form of coarse-graining procedure must be carried out, which essentially 
amounts to treating two mechanical states as the same microstate if the differences in their 
position and momentum variables are 'small enough'. The partition function then takes the 
form of an integral. For instance, the partition function of a gas consisting of A'^ molecules 
is proportional to the 6A/'— dimensional phase-space integral. 



Z(/3)~/ d^pid^x'exp[-pH{pi,x% 



where H = H{pi,x'), {i = 1, ...,N) is the classical Hamiltonian (total energy) function. 

More generally, the so-called configuration integral, as used in probability theory, infor- 
mation science and dynamical systems, is an abstraction of the above definition of a partition 
function in statistical mechanics. It is a special case of a normalizing constant in probability 
theory, for the Boltzmann distribution. The partition function occurs in many problems of 
probability theory because, in situations where there is a natural symmetry, its associated 
probability measure, the Gibbs measure (see below), which generalizes the notion of the 
canonical ensemble, has the Markov property. 

Given a set of random variables Xi taking on values , and purely potential Hamiltonian 
function H{x^), {i = 1, N), the partition function is defined as 

Z(/3) = ^exp[-/3i?(x^)]. (55) 

The function H is understood to be a real- valued function on the space of states {Xi , X2 , • • • } 
while /? is a real- valued free parameter (conventionally, the inverse temperature). The sum 
over the x^ is understood to be a sum over all possible values that the random variable Xi 
may take. Thus, the sum is to be replaced by an integral when the Xi are continuous, rather 
than discrete. Thus, one writes 



Z{l3) = jdx'exp [-PH{x')] 



for the case of continuously- varying random variables Xi . 

The Gibbs measure of a random variable Xi having the value is defined as the proba- 
bility density function 

P(X. =x') = ^ exp \-0E(x')] = [-/^^(^^)] 
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where E{x^) — H{x^) is the energy of the configuration . This probabihty, which is now 
properly normahzed so that < P{x^) < 1, can be interpreted as a hkehhood that a specific 
configuration of values a;*, (z = 1,2, ...N) occurs in the system. P(x') is also closely related 
to f2, the probability of a random partial recursive function halting. 

As such, the partition function Z(f3) can be understood to provide the Gibbs measure on 
the space of states, which is the unique statistical distribution that maximizes the entropy 
for a fixed expectation value of the energy, 

= — dp — • 

The associated entropy is given by 



representing 'ignorance' + 'randomness'. 

The principle of maximum entropy related to the expectation value of the energy {H) , is 
a postulate about a universal feature of any probability assignment on a given set of propo- 
sitions (events, hypotheses, indices, etc.). Let some testable information about a probability 
distribution function be given. Consider the set of all trial probability distributions which 
encode this information. Then the probability distribution which maximizes the informa- 
tion entropy is the true probability distribution, with respect to the testable information 
prescribed. 

Applied to the crowd dynamics, the Boltzman's theorem of equipartition of energy states 
that the expectation value of the energy {H) is uniformly spread among all degrees-of- freedom 
of the crowd (that is, across the whole crowd manifold M). 



6.3 Free energy, Landau's phase transitions and Haken's synergetics 

All thermodynamic-like properties of a multi-component system like a human (or robot) 
crowd may be expressed in terms of its free energy potential, T = —ksThiZlP), and its 
partial derivatives. In particular, the physical entropy S of the crowd is defined as the 
(negative) first partial derivative of the free energy J- with respect to the control parameter 
temperature T, i.e., S = —OtT, while the specific heat capacity C is the second derivative, 
C = TdrS. 

A phase of the crowd denotes a set of its states that have relatively uniform psycho- 
physical properties. A crowd phase transition represents the its transformation from one 
phase to another (see e.g., [46ll56j ). In general, the crowd phase transitions are divided into 
two categories: 

• The first-order phase transitions, or, discontinuous phase transitions, are those that 
involve a latent heat C . During such a transition, a crowd either absorbs or releases a 
fixed (and typically large) amount of energy. Because energy cannot be instantaneously 
transferred between the system and its environment, first-order crowd transitions are 
associated with mixed-phase regimes in which some parts of the crowd have completed 
the transition and others have not. This forms a turbulent spatioi-temporal chaotic 
interphase, difficult to study, because its dynamics can be violent and hard to control. 
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• The second-order phase transitions are the continuous phase transitions, in the entropy 
5* is continuous, without any latent heat C . They are purely entropic crowd transitions, 
which are at the focus of the present study. 

In Landau's theory od phase transitions (see [ISHSS]), the probability density function P 
is exponentially related to the free energy potential J^, i.e., P ~ e~-^^'^\ if is considered as 
a function of some order parameter o. Thus, the most probable order parameter is determined 
by the requirement = min. Therefore, the most natural order parameter for the crowd 
dynamics would be its entropy S. 

The following table gives the analogy between various systems in thermal equilibrium 
and the corresponding nonequilibrium systems analyzed in Haken's synergetics |24[ 1251 126) : 



System in thermal equilibrium 


Nonequilibrium system 


Free energy potential J- 


Generalized potential V 


Order parameters Oi 


Order parameters oi 


aoi 


a. = 

OOi 


Temperature T 


Control input u 


Entropy S 


System output y 


Specific Heat c 


System efficiency e 



In particular, in case of human biodynamics [461 156j. natural control inputs Ui are mus- 
cular forces and torques, Fi, natural system outputs yi are joint coordinates and mo- 
menta Pi, while the system efficiencies represent the changes of coordinates and mo- 
menta with changes of corresponding muscular torques for the ith active human joint, 

q _ dq^ p _ dpi_ 

dFi ' dFi ■ 

6.4 Heat equation, Dirichlet action and gradient flow on a Rieman- 
nian manifold 

The heat equation 

u = Au, (56) 

on a compact Riemannian manifold M with static metric {dtg = 0), where u : [0, T] x Af — ^ R 
is a scalar field, can be interpreted as the gradient flow for the Dirichlet action 

Eiu) ■.= \ [ |Vu|2 d/x, (57) 

using the inner product, {ui,U2)^i ■= Jjyj itiU2 d^, associated to the volume measure dfj,. This 
can be proved if we evolve u in time at some arbitrary rate it, an application of integration 
by parts formula, 

[ uX/^X" dn=- [ (Vau)X" dfi 

JM JM 

(where div(X) := V^X" is the divergence of the vector-field X", which validates the Stokes 
theorem, /j^^div(X) d^ — 0), gives 

dtE{u)^~ {Au)u d^i ^ {~Au,u)f,, (58) 

JM 
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from which we see that (|56l) is indeed the gradient flow for ((58| with respect to the inner 
product. In particular, if u solves the heat equation (j56p . we see that the Dirichlet energy is 
decreasing in time, 

dtE{u) = - /" |Aup dfi. (59) 
Jm 

Thus we see that by representing the parabolic PDE (|56p as a gradient flow, we automatically 
gain a controlled quantity of the evolution, namely the energy functional that is generating 
the gradient flow. This representation also strongly suggests that solutions of (|56p should 
eventually converge to stationary points of the Dirichlet energy (|57p . which by (j58p are 
harmonic functions (i.e., the functions u with Au = 0). As an application of the gradient 
flow interpretation, we can assert that the only periodic (or, "breather") solutions to the 
heat equation (j56p are the harmonic functions (which must be constant if the manifold M 
is compact). Indeed, if a solution u was periodic, then the monotone functional E must be 
constant, which by (j59p implies that u is harmonic as claimed. 



6.5 Lyapunov exponents and Kolmogorov— Sinai entropy 

A branch of nonlinear dynamics has been developed with the aim of formalizing and quan- 
titatively characterizing the general sensitivity to initial conditions. The largest Lyapunov 
exponent A, together with the related Kaplan-Yorke dimension dxY and the Kolmogorov- 
Sinai entropy hx s are the three indicators for measuring the rate of error growth produced 
by a dynamical system [121 SHI \EE\ ■ 

The characteristic Lyapunov exponents are somehow an extension of the linear stability 
analysis to the case of aperiodic motions. Roughly speaking, they measure the typical rate 
of exponential divergence of nearby trajectories. In this sense they give information on the 
rate of growth of a very small error on the initial state of a system [HI |9] . 

Consider an nD dynamical system given by the set of ODEs of the form 

i - fix), (60) 

where x = (xi,...,a;„) G R" and / : M" — * R". Recall that since the r.h.s of equation 
((60|l does not depend on t explicitly, the system is called autonomous. We assume that / is 
smooth enough that the evolution is well-defined for time intervals of arbitrary extension, 
and that the motion occurs in a bounded region R of the system phase space M. We intend 
to study the separation between two trajectories in M, x{t) and x'{t), starting from two 
close initial conditions, x{0) and a;'(0) = x{0) + Sx{0) in C M, respectively. 

As long as the difference between the trajectories, 6x{t) = x'{t) —x{t), remains infinitesi- 
mal, it can be regarded as a vector, z{t), in the tangent space T^M of M. The time evolution 
of z(t) is given by the linearized differential equations: 



z,{t). 

x{t) 



Under rather general hypothesis, Oseledets [SS] proved that for almost all initial conditions 
x{0) S R, there exists an orthonormal basis {e^} in the tangent space T^M such that, for 
large times, 

z{t) = CiBi exp(Aii), (61) 
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where the coefficients {ci} depend on ^(O). The exponents Ai > A2 > ■ • • > are called 
characteristic Lyapunov exponents. If the dynamical system has an ergodic invariant measure 
on M, the spectrum of LEs {Xi} does not depend on the initial conditions, except for a set 
of measure zero with respect to the natural invariant measure. 

Equation (pT|) describes how an nD spherical region i? = 5" C M, with radius e centered 
in a;(0), deforms, with time, into an ellipsoid of semi-axes ei(t) — eexp(Ait), directed along 
the Ci vectors. Furthermore, for a generic small perturbation 6x{0), the distance between 
the reference and the perturbed trajectory behaves as 

\dx{t)\ - |,5a;(0)| exp(Ait) [1 + O (exp -(Ai - As)*)] . 

If Ai > we have a rapid (exponential) amplification of an error on the initial condition. In 
such a case, the system is chaotic and, unpredictable on the long times. Indeed, if the initial 
error amounts to Sq = \Sx(0)\, and we purpose to predict the states of the system with a 
certain tolerance A, then the prediction is reliable just up to a predictability time given by 

This equation shows that Tp is basically determined by the positive leading Lyapunov expo- 
nent, since its dependence on 60 and A is logarithmically weak. Because of its preeminent 
role, Ai is often referred as 'the leading positive Lyapunov exponent', and denoted by A. 

Therefore, Lyapunov exponents are average rates of expansion or contraction along the 
principal axes. For the ith principal axis, the corresponding Lyapunov exponent is defined 
as 

A, = lim {{1/t) ln[L,(i)/L,(0)]}, (62) 

t — ^00 

where Li{t) is the radius of the ellipsoid along the ith principal axis at time t. 
An initial volume Vq of the phase-space region Rq evolves on average as 

V{t) = yoe(^i+^2+-+A.„)t^ (g3) 
and therefore the rate of change of V{t) is simply 



T/(i) = ^A.y(t). 



i=l 



In the case of a 2D phase area A, evolving as A{t) = Aoe'^^^'^^^^*' , a Lyapunov dimension 
is defined as 

-d(ln(7V(e)))- 



= lim , , 

e^o [ d(ln(l/e)) 

where iV(e) is the number of squares with sides of length e required to cover A{t), and d 
represents an ordinary capacity dimension, 

dc — lim \ T—-—— 

e-.o Vln(l/e) 
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Lyapunov dimension can be extended to the case of nD phase-space by means of the 
Kaplan-Yorke dimension [621 1851 as 

, Ai + A2 + • • • + Aj 

dxY =J + TT . 

|Aj+i| 

where the Ai are ordered (Ai being the largest) and j is the index of the smallest nonnegative 
Lyapunov exponent. 

On the other hand, a state, initially determined with an error (5x(0), after a time enough 
larger than 1/A, may be found almost everywhere in the region of motion R e M. In this 
respect, the Kolmogorov-Sinai (KS) entropy^ h^Sj supplies a more refined information. The 
error on the initial state is due to the maximal resolution we use for observing the system. 
For simplicity, let us assume the same resolution e for each degree of freedom. We build 
a partition of the phase space M with cells of volume e"*, so that the state of the system 
at i = <o is found in a region Rq of volume Vq = e'^ around x{to). Now we consider the 
trajectories starting from Vq at and sampled at discrete times tj = j t {j = 1,2,3, ... ,t). 
Since we are considering motions that evolve in a bounded region R C M, all the trajectories 
visit a finite number of different cells, each one identified by a symbol. In this way a unique 
sequence of symbols {s(0), s(l), s(2), . . . } is associated with a given trajectory x{t). In a 
chaotic system, although each evolution x(t) is univocally determined by a; (to), a great 
number of different symbolic sequences originates by the same initial cell, because of the 
divergence of nearby trajectories. The total number of the admissible symbolic sequences, 
N(e, t), increases exponentially with a rate given by the topological entropy 

hr = lim lim — IniVfe, t) . 

e^Oi^oo t 

However, if we consider only the number of sequences Neff{e, t) < N{e, t) which appear 
with very high probability in the long time limit - those that can be numerically or exper- 
imentally detected and that are associated with the natural measure - we arrive at a more 
physical quantity called the Kolmogorov-Sinai (or metric) entropy, which is the key entropy 
notion in crgodic theory [16] : 

hxs = lim lim i InTVe// (e, i) < /it- (64) 

hKS quantifies the long time exponential rate of growth of the number of the effective coarse- 
grained trajectories of a system. This suggests a link with information theory where the 
Shannon entropy measures the mean asymptotic growth of the number of the typical se- 
quences - the ensemble of which has probability almost one - emitted by a source. 

We may wonder what is the number of cells where, at a time t > to, the points that 
evolved from Rq can be found, i.e., we wish to know how big is the coarse-grained volume 
V{e, t), occupied by the states evolved from the volume Vq of the region Rq, if the minimum 
volume we can observe is Vmin — f^. As stated above we have 



T/(t) -Foexp(t^A,). 



i=l 
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However, this is true only in the hmit e ^ 0. In this (unreahstic) hmit, V{t) = Vq for 
a conservative system (where — ^) ^(^) < ^ ^o^' ^ dissipative system (where 

X]f=i -^i < 0)- As a consequence of hmited resohition power, in the evolution of the volume 
Vq — e'^ the effect of the contracting directions (associated with the negative Lyapunov 
exponents) is completely lost. We can experience only the effect of the expanding directions, 
associated with the positive Lyapunov exponents. As a consequence, in the typical case, the 
coarse grained volume behaves as 

y(e,i) ~ Vbe^^^.>»^'^*, 

when Vo is small enough. Since Neff{e,t) oc F(e,i)/Vo, one has: hxs = J2x >o^i- "^^^^ 
argument can be made more rigorous with a proper mathematical definition of the metric 
entropy. In this case one derives the Pesin relation [751 [TB]: Iiks < X^a >o Because of 
its relation with the Lyapunov exponents, or by the definition (|64p , it is clear that also hxs 
is a fine-grained and global characterization of a dynamical system. 

The metric entropy is an invariant characteristic quantity of a dynamical system, i.e., 
given two systems with invariant measures, their KS-entropies exist and they are equal iff 
the systems are isomorphic [5]. 

Finally, the topological entropy on the manifold M equals the supremum of the Kolmogorov- 
Sinai entropies, 

h{u) = sup{/iKs(w) = hf^{u) : ^ e P„(A/)}, 

where u : M — > M is a continuous map on M, and fj, ranges over all u— invariant (Borel) 
probability measures on M. Dynamical systems of positive topological entropy are often 
considered topologically chaotic. 
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